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Abstract 

Background: Human Rhinoviruses (HRVs) have high genetic diversity and three species have been described: 
HRV-A, HRV-B, and the recently recognized HRV-C, which has been rapidly identified worldwide. 

Findings: In the present study, we report the frequency and diversity of Human Rhinovirus (HRV) strains 
circulating in Panama from children hospitalized with respiratory infections. 

Conclusions: HRVs of species A, B and C have been identified with a predominance of HRV-A and HRV-C over 
HRV-B, and marked genetic diversity within each species. 
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Findings 

Human rhinoviruses (HRVs) are the most common 
causative agents of upper respiratory tract infections, 
but are also associated with more severe diseases such as 
pneumonia or acute wheezing related to bronchiolitis 
and acute asthma in children [1-4]. HRV infection 
occurs in all age groups and is responsible for 25% to 
50% of respiratory infections presenting as influenza-like 
illnesses [5-7]. 

HRVs have been classified into the genus Enterovirus, 
family Picornaviridae [8]. As other picornaviruses, they 
are small, non-enveloped viruses with a 7200 bp single- 
stranded, positive-sense mRNA genome and a long open 
reading frame encoding four capsid proteins VP4, VP2, 
VP3 and VP1, and seven non-structural proteins 2A, 2B, 
2C, 3A, 3B, 3C and 3D. 

HRVs have high genetic diversity and three species 
have been described: HRV-A, HRV-B, and the recently 
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recognized HRV-C [9], which has been rapidly identified 
worldwide [10-12]. 

Reports that infection with different HRV species re- 
sult in different clinical outcomes are controversial: 
some show correlations between a given viral serotype 
or species and its capacity to invade the upper or 
lower respiratory tract [13], whereas others provide 
evidence for a more frequent role for HRV-C in lower 
respiratory tract infections, associated with more severe 
disease [9,11]. In a recent study, HRV rates were high 
among hospitalized children and the elderly, but HRV 
was also detected among asymptomatic children. HRV- 
A and HRV-C were associated with illness requiring 
hospitalization [14]. 

Very little in general is known about the genetic diver- 
sity of respiratory viruses in Central America or Panama, 
nor specifically about rhinoviruses. In this study we inves- 
tigate the genetic variability of Rhinovirus from children 
under five years, hospitalized in Panama with respiratory 
infections. Stored nasal swab samples were used in this 
study and the protocol was approved by Ethics Committee 
of Gorgas Institute. Sixty-two nasal swab samples positive 
for HRV by real time PCR [15] from a total 118 samples 
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from 420 study subjects, collected from August 2010 («; 
178) to June 2011 (n: 242), were analyzed by RT-PCR and 
sequencing. Briefly, the viral RNA was extracted using 
QIAamp Viral RNA Mini Kit (Qiagen, GmbH, Hilden, 
Germany) according to the manufacturer's instructions. 
Amplification of a 542 bp fragment containing the HRV 
VP4 and partial VP2 region was performed by RT-PCR 
using a reverse primer identical to that described previ- 
ously [16] and external and semi-nested forward primers 



designed by Xiaoyan Lu and Dean Erdman at CDC (un- 
published data). cDNA synthesis and PCR were carried 
out using QIAGEN OneStep RT-PCR Kit and Taq PCR 
Master Mix Kit and PCR products were purified using 
a Wizard® SV Gel and PCR Clean-Up System (Promega 
Corporation Madison, WI 53711 USA). Sequencing was 
performed in both directions using the same primers 
as in the PCR with an ABI Prism BigDye Terminator 
Cycle Sequencing Ready Reaction Kit according to the 
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Figure 1 Phylogenetic analysis of the VP4/VP2 gene from the 62 HRVs of this and reference strains for each species. Phylogenetic trees 
were constructed by maximum likelihood by using PhyML nucleotide substitution model was HKY + I + G, aLRT supports are depicted above the 
nodes. Virus sequence names shown symbols were reference strains of each HRV species with GenBank accession numbers in parentheses. The 
GenBank accession numbers assigned for the Panamanian samples were JN797533 to JN797594. 



Franco et al. Virology Journal 201 2, 9:257 
http://www.virologyj.eom/content/9/1/257 



Page 3 of 6 



manufacturer's instructions on a 3130XL Genetic 
Analyzer (Applied Biosystems, Foster City, CA). 

Alignment and sequence analyses were performed 
using Bioedit, ClustalW and MEGA5 software. The most 
suitable model for nucleotide substitution was estimated 
with Modelgenerator. Phylogenetic trees were con- 
structed under maximum likelihood criteria and branch 
supports were calculated by the approximate likelihood 
ratio test (aLRT). One hundred and fifty published HRV 
sequences within the VP4-VP2 region were obtained 
from the GenBank database (NCBI). The VP4/VP2 re- 
gion was used to analyze the sequence variability within 
the much larger data set of published sequences of all 
HRV-A, HRV-B serotypes, and also from the recently 
described species C variants, which have been assigned 
to 28 provisional types [12]. 

Alignment of the sequences in the present study with 
the HRV reference strains resulted in the three HRV-A, 
HRV-B, and HRV-C genetic clusters (Figure 1). Higher 
frequencies of detection were observed for HRV-A (60%, 
n = 37) and HRV-C (32%, n = 20), over HRV-B species 
(8%, n = 5). In a recently published study of both chil- 
dren and adults HRV-A was reported as more common 
than HRV-C in children aged less than 1 year, while 
HRV-C was the most common species in juvenile 
patients aged 1 to 19 years [14]. In our study, HRV-A 
was the most common species in children aged less than 
1 year (24 of 44 samples), but also in patients aged 1 to 
5 years (13 of 18 samples). 

The same region of VP4-VP2 region for each positive 
sample amplified and sequenced was used to identify the 
serotypes within each species. To do this, we constructed 
phylogenetic trees with the dataset of published sequences 
belonging to all HRV-A, HRV-B serotypes, and also from 
the recently described HRV-C species variants [9,12]. 

Variants in HRV-C species are not currently assigned 
to serotypes. However, some attempts to assign types 
with temporary names based on tree position have been 
performed [11], and more recently a proposal for the 
designation of HRV-C provisional types 1-28 has been 
published, based only on VP4/VP2 sequences [12]. Fol- 
lowing the latter classification, the 28 published 
sequences were used to construct a phylogenetic tree 
with the 20 HRV-C samples identified in this work. 
Fourteen sequences were assigned to previously reported 
types, however, six samples (308526, 308759, 308769, 
308595, 308623 and 308248) could not be assigned to 
any of the 28 provisional types (Table 1, Figure 2-C). 
These samples exhibited a percentage of p-distance 
greater than 10% with respect to the 28 HRV-C types. 
Thus, following the criterion proposed by Simmonds 
et al., (> 10% of divergence in VP1-VP4 region between 
different types) the untyped strains could be assigned to 
new types. 



Table 1 HRV serotypes detected among study subjects 
within each species 

Species and serotype No. of samples with the indicated serotype 
Species A 



A12 4 

A57 1 

A24 5 

A71 4 

A58 1 

A9 1 

A40 4 

A96 2 

NS1 4 

NS2 3 

A60 1 

A32 2 

A77 1 

A59 1 

A20 1 



Untyped (308639) 1 

Untyped (308868) 1 
Species B 

B3 1 

B6 1 

B52 2 

B97 1 
Species C 

Pat8 2 

Patl 1 1 

Pat24 1 

Pat9 1 

Pat7 2 

Pat22 1 

Patl 8 1 

Pat3 1 

Pat19 2 

Pat25 1 

Patl 5 1 

Untyped (308526) 1 

Untyped (308759) 1 

Untyped (308769) 1 

Untyped (308595) 1 

Untyped (308623) 1 

Untyped (308248) 1 

Total 62 
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Figure 2 (See legend on next page.) 
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(See figure on previous page.) 

Figure 2 Phylogeny of the VP4-VP2 sequences amplified of HRV variants detected in the study in comparison of those sequences from 
serotype reference strains of species A and B and provisionally assigned for species C. Phylogenies were constructed under maximum 
likelihood criterion (PhyML software), with the following models of nucleotide substitution: GTR + I + G for species A and B and HKY+ I + G for 
species C. The GenBank accession numbers assigned for the Panamanian samples were JN797533 to JN797594. For species B and C variants, 
HRV-A (GenBank accession number DQ473509) was used to root the trees; while sequence HRV-B14 (GenBank accession number NC_001490) 
was used to root the species A tree. 



Table 1 and Figure 2- A, show that of 37 HRV-A sam- 
ples, 35 were assigned to previously reported serotypes, 
including 7 samples clustering with the recently assigned 
new HRV-A serotypes, named NS1 and NS2 [11]. Two 
samples (308639 and 308868) could not be properly 
assigned to the actually defined serotypes. Sequence 
308639 displayed a p-distance greater than 10% with re- 
spect to all of the HRV-A serotypes, except for serotypes 
29 and 44. Because of this, and following the criterion 
postulated by Simmonds et al., for HRV-C (> 10% of di- 
vergence between types) [12], this strain should not be 
assigned to a new serotype. Sequence 308868 shows a p- 
distance greater than 10% with respect to every other 
serotype. Thus, following the same criterion, this strain 
could be provisionally assigned to a new serotype. 

Five samples identified as HRV-B were clustered into 
four previously reported HRV-B serotypes (Table 1, 
Figure 2). 

The HRV strains analyzed showed substantial genetic 
diversity, a fact that is supported by a Shannon entropy 
analysis performed with the sequences from Figure 1, 
which displayed values between 0.6 and 1 for the major- 
ity of the sites analyzed (data not shown, available on 
request). 

Marked genetic diversity within HRV-A and HRV-C 
over HRV-B species has been demonstrated in this 
study, with HRV-B sequences being the most homoge- 
neous species, showing an intra-species nucleotide p- 
distance of 0.19, followed by rhinovirus A, with 0.20. 
The most heterogeneous was species C, with a p- 
distance of 0.25. 

The emergence of previously untyped serotypes of 
HRVs in our study represents a challenge for HRV 
researches to develop the necessary tools to allow for re- 
liable classifications. 

Additional studies involving subject of all ages hospita- 
lized and also ambulatory will be done in Panama, in 
order to elucidate the relationships between HRV species 
and to define the epidemiological profile and genetic 
characteristics of each one. 

Panamanian strains were deposited in the GenBank 
Database under accession numbers [Gen-Bank: JN797533 
to JN797594]. 
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